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ABSTRACT 

M82 X-I is the brightest ultraluminous X-ray source in starburst galaxy M82 and is one of the best 
intermediate mass black hole candidates. Previous studies based on the Rossi X-ray Timing Explorer 
/Proportional Counter Array (RXTE/PCA) reported a regular X-ray flux modulation of M82 with 
a period of 62 days, and attributed this periodic modulation to M82 X-1. However, this modulation 
is not necessarily from M82 X-1 because RXTE/PCA has a very poor spatial resolution of 1°. 
In this work, we analyzed 1000 days of monitoring data of M82 X-1 from the Swift/X-iay telescope 
(XRT), which has a much better spatial resolution than RXTE/PCA. The periodicity distribution 
map of M82 reveals that the 62-day periodicity is most likely not from M82 X-1, but from the summed 
contributions of several periodic X-ray sources 4" southeast of M82 X-1. However, Swift/XPX is not 
able to resolve those periodic sources and locate the precise origin of the periodicity of M82. Thus, 
more long-term observations with higher spatial resolution are required. 

Subject headings: galaxies: individual (M82) — X-rays: binaries 


1. INTRODUCTION 

Ultraluminous X-ray sources (ULXs; Fabbiano 2006) 
discovered in other galaxies exhibit high X-ray luminosi¬ 
ties Lx > 2 X 10^® erg and are possibly the long 
sought after intermediate mass black holes (IMBHs) radi¬ 
ating at sub-Eddinton levels. Recent studies have shown 
that, except for a few cases, many ULXs can be stellar 
mass black holes i n an ultraluminous state un der super¬ 
critical accretion (I Gladstone fc Robertsll200^ . Indeed, 
optical monitori ng campaigns of two ULXs, MlOl ULX- 1 
(ITT] et a1 .11201 .in and NGC 7793 P13 (IMotch et aHl^oTl . 
have confirmed both to be stellar mass black holes, while 
NuSTAR observations of M82 have discovered that a 
transient ULX, M82 X-2, is actually pow ered by an ac¬ 
creting neutron star (|Bachetti et al.l[2014D . 

M82 X-1 (CXO J095550.2-h694047) is the brightest 
ULX in the starburst galaxy, M82, at a distance of 
3.6Mpc (iFreedman et all 1199^ . and arguably the best 
IMBH candidate in the local universe. M82 X-1 ex¬ 
hibits extremely high luminosities t hat can reach 10^^^ 
erg in the 0.5 — 10 keV bands (iKaaret et ahll^OOll : 
iMatsumoto et al.l 120011 : IKaaret et al.l 1200911 . suggesting 
an IMBH of about 10^ Mq if radiating at the Eddington 
level. Recently, Pasham et al. (2014) re-analyzed 6- 
year Rossi X-ray Timing Explorer {RXTE) X-ray obser¬ 
vations and revealed the high-frequency, 3:2 ratio, twin- 
peak quasi-periodic oscillations (QPOs) of 3.3 Hz and 
5Hz, which, in combination with the low-frequency QPO 
revealed by an XMM-Newton observation, suggests an 
IMBH of 415 ± 63M(^ un der the relativistic precession 
model (|Motta et al.ll20l4l . 

Interestingly, X-ray observations of M82 with the Pro- 
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portional Counter Array (PCA) on board RXTE over 
eight months has revealed regular m odulations of about 
20% with a per i od of about 62 days (IKaaret et al.ll2006l: 
IKaaret fc Fend I2007L hereafter KF07). IKaaret et all 
( 200611 interpreted such a period as the orbital period 
for M82 X-1, and consequently derived a giant or su¬ 
pergiant companion for M82 X-1. However, the analy¬ 
sis of another 4.5-year RXTE/PCA light curve of M82 
has shown a sudden phase shift of the 62-day modula¬ 
tion, suggesting the 62-day X-ray period may be caused 
by a precessing accretion di sk around the black hole 
(IPasham &: Strohmaver|[2M^ hereafter PS13). 

While previous studies have attributed the periodic 
modulations in the RXTE/PCA light curves to M82 X-1, 
the brightest X-ray source in M82, it is not necessarily 
true because RXTE/PCA has a very poor spatial resolu¬ 
tion of ~ 1°, and the flux includes contributions from all 
X-ray sources in M82 (11' x 4') and nearby galaxies. In 
contrast, the X-ray telescope (XRT) on board the Swift 
satellite has a much better spatial resolution with a half 
power diameter of 18", and may resolve the bright point 
sources in M82. In this Letter, we use Swift/XRT obser¬ 
vations over 1000 days to scrutinize the X-ray periodicity 
of M82. The data analysis and results are presented in 
§2, and the discussion follows in §3. 

2. DATA ANALYSIS AND RESULTS 

The Swift/XRT observations of M82 operated in the 
photon-counting mode are used in this work, which are 
retrieved from the HEASARC archive. These include 
a total of 106 observations spanning from 2012 April 
to 2015 January. 70 observations in 2014 are not used 
because they are badly contaminated by the supernova 
SN 2014J in M82. Inspection of these XRT observations 
confirms that X-ray photons from M82 X-1 and nearby 
sources dominate the field within ^ 30' of M82. 

2.1. Stacked images of M82 with Swift observations 

There are three other ULXs near M82 X-1, 

e.g., M82 X-2, CXOU J095551.2 -L694044 (here- 
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after X - 3 ) and CXOU J 09555 0 6 +694944 (here¬ 
after X-4; L iu fc Mirabell 120051 : iKong et all 120071 : 
IStrohmaver fc Mushotzk'vll2003n \^th 2" — 6" separations 
from each other. Those four ULXs can only be clearly 
resolved by the Chandra X-ray observatory (see Figure 
1(a)), unfortunately there are not enough Chandra obser¬ 
vations of M82 X-1 in order to perform periodicity analy¬ 
sis. In contrast, Swift/XRT has more than a hundred ob¬ 
servations of M82 X-1, and it can marginally resolve the 
four ULXs in a single observation (Figure 1(b)), though 
the claimed half power diameter of the point spread func¬ 
tion of XRT is 18". The image spatial resolution of XRT 
is one pixel (2".36), but the positional accuracy from cen- 
troiding can be higher when there are ample photons. To 
reveal of all the faint and bright X-ray sources with more 
accurate positions, we stacked the total 106 XRT obser¬ 
vations. Before image stacking, one pixel was equally 
divided into 10 x 10 parts, and each part (a sub-pixel) 
was assigned 1% of the pixel counts. Then the total 106 
images in sub-pixels were matched under WCS coordi¬ 
nates and stacked into one (Figure 1(c)). 

This stacked image (Figure 1(c)) reveals that the four 
ULXs are not blurred into one unresolved source, but an 
asymmetric structure in which X-1 and the other three 
ULXs are separated. X-2 is the second brightest X-ray 
source in M82, but it is not evident in Figure 1(c) because 
it is a transient and, for most of the time, it is in a qui¬ 
escent state. Among the 106 Swift observations, only a 
dozen observations in which X-2 is brighter than X-3 and 
X-4. X-3 is also a variable X-ray source, and is visible in 
46 observations . X-4 was reported as a supernova rem - 
nant candidate (|Kong et al.ll2007HChiang fc Kong|l2011h . 
which is supposed to have relatively constant flux. 

Furthermore, to inspect whether XRT can resolve X-1 
from X-3 or not, 46 observations in which X-3 is brighter 
than the nearby background were selected , and stacked 
into one image (Figure 1(d)). The rest of the 60 observa¬ 
tions are stacked and shown in Figure 1(e). These figures 
imply that X-1 can be roughly resolved from X-3, while 
counts from X-2, X-3, and X-4 within 4" mix up with 
each other and can hardly be resolved by XRT. Thus, 
sources with separations larger than 4" can marginally 
resolved by XRT, but not vice versa. 


2.2. Phase offset of the 62 day X-ray Period 

Based on the long-term observations, Swift/XKT data 
can be used to test the 62-day modulation (KF07) of 
M82. Source counts were extracted from a circular re¬ 
gion 18" in ra dius centered on X- 1, with 73% of en¬ 
ergy encircled (|Moretti et al.l [200^ . The background 
was defined by a circular region with a radius of 20 pix¬ 
els (47") placed near M82 X-1, but avoiding the con¬ 
tamination from other point sources. The 0.3-10 keV 
background-subtracted light curve, binned per observa¬ 
tion, was used to look for the periodicity of X-1 with the 
Lomb-Scargle(L-S) algorithm (iLomblllQT^ iScarglel[19821 : 
?). The power of the periodo gram was normalized b y 
the total variance of the data (|Horne fc Baliun^ll986ll . 
As shown in Figure [5] (a), the main peak appears at the 
period of 61.8 ± 2.9 days, and the error is estimated by 
calculatin g the F'WHM of the main peak in the peri¬ 
odogram (|Kovacslll981ll . 

The significance of the main peak in periodogram can 
be given by the false alarm probability (hereafter fap) as¬ 


suming all of the datapoints were generated by stochastic 
noises. The fap based on the Gaussian white noise of the 
main peak is less than 0.001 (dashed horizontal lines in 
Figure 2(a)), i.e., the confidence level of this period gen¬ 
erated by a signal instead of the white noise background 
is more than 99.9%. While significant peaks under the 
white noise background may have low significance when 
taking the red noise component into consideration. In 
contrast to white noise, whose power is independent of 
frequency, red noise exhibits a power-law spectrum in 
the form of P{ijj) oc (/3 = 0 for white noise), which 
can lead to relatively large a periodic peaks in th e power 
spectrum at low frequencies (iGruber et al.ll20lHl . 

In actuality, the periodogram (Figure 2(a)) does not 
show evidence of increasing slope for red noise at long 
periods, but is rather flat wit h a best-fit power law in¬ 
dex of —0.32 ([Vaughanl 1200511 . In spite of this, we es¬ 
timated the significance of the observed data under a 
red noise background. 100,000 red noise synthetic light 
curves with a mean and varianc e equal to the ones of 
the actual data were ge nerated ([Davies fc Harte 1 119871 : 
iTimmer fc Koenig|ll995f) . Each synthetic light curve has 
the same sampling time series as the actual data set and 
was u sed to conduct L-S periodograms (see lKaaret et al.l 
120061 for more details). We found only two cases that 
power values at a period of 61.8 ± 2.9 days excess the 
power (12.24) of the observed highest peak. Namely, the 
probability that a red noise produces a peak at the given 
period ranges higher than 12.24 is only 2 x 10“®. The fap 
of the highest peak at the inspected frequencies (40-80 
days) based on red n oise is 0.0012 (confidence of 99.88%) 
(iGharisi et al.ll2015il . Hence, the 62-day periodicity from 
Swift/XKT data is significant enough to be produced by 
a real periodic signal, and confirms the ~ 62—day mod¬ 
ulation of M82 discovered by KF07 with RXTE/PCA. 

The phase shift of the 62-day period was reported by 
PS13 utilizing RXTE/PCA observations. Following the 
work in PS13, we folded the new Swift/XPCP light curve 
(MJD 56022.9-MJD 57053.4) at 62 days to test the sta¬ 
bility of the ^ 62—day periodicity. The previous RXTE 
light curves from PS 13 consist of two segments, which 
cover MJD 53250.6-54225.4 (segment 1, the same data 
used in KF07) and MJD 54353.8-55195.8 (segment 2), 
respectively. Start times of the Swift light curve and the 
two RXTE light curves were set to be the same (MJD 
53250.6). Then the three light curves were folded at the 
period of 62 days (see Figure IHb)). The two RXTE 
light curves indeed reveal a 0.4 phase offset as reported 
in PS13. Compared with RXTE light curves, the Swift 
light curve shows changes of about 0.2 in phase relative 
to RXTE segment 1, and 0.4 in phase relative to RXTE 
segment 2. It confirms the result in PS13 that the 62-day 
X-ray periodicity of M82 has phase changes, indicating 
that it is not stable. 

2.3. Origin of the periodicities in M82 

Since the 18" radius circular region used in §2.1 en¬ 
circles not only X-1 but many other X-ray sources, it is 
not clear whether the 62-day periodicity originates from a 
certain single source, like X-1, or from summed contribu¬ 
tions of several X-ray sources. Hence, we used a smaller 
source region with a radius of 4", which can separate X-1 
from other X-ray sources, yet contains 30% of the energy 
encircled in the 18" region used in §2.2, to construct the 
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Fig. 1.— X-ray images of M82 from Chandra/AQIS and 5m/t/XRT. Each panel is shown on the same spatial scale, and north is up. The 
spatial size of (a), (b), (d), and (e) is 18^' x 18^^ and the size of (c) is 36'' x 36". The positions of X-1, X-2, X-3, and X-4 are marked by 
crosses and corresponding numbers, (a) Chandra/KQ\S image from observation ObsID 5644. (b) Swift/X.'KT image from observation ObsID 
00032503069, which shows that the four ULXs can be marginally resolved by Swift/'KB^. (c) The stacked image of all 106 5'iu2/i/XRT 
observations used in this study. The counts concentrated on X-1 and X-4 and reveal an asymmetric structure, (d) The stacked image of 46 
Swift/X.'KT observations in which X-3 is obviously brighter than the background, (e) The stacked image of the remaining 60 observations. 
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Fig. 2.— (a) L-S periodogram of a Swift/'K.'KT circular region centered at X-1 with a radius of 18". A peak appears at 61.8 di 2.9 days. 
Dashed, horizontal lines denote the fap of 0.5, 0.1, 0.05, 0.01, 0.005, and 0.001 (or confidence levels of 50%, 90%, 95%, 99%, 99.5%, and 
99.9%), based on the Gaussian noise background, (b) Normalized phase light curves of Swift^ RXTE segment 1, and RXTE segment 2 
folded at 62 days. Each point is the averaged count rate of the observations falling within the given phase bin, and the error bar is shown 
as the standard deviation divided by the square root of the number of points in that bin. The solid sinusoidal curves are the best fourier 
fit of the data. 


periodicity map of M82 and analyze the origin of the 62- 
day modulation. We divided the region around M82 X-1 
into a grid with steps of 0.5 pixel (i.e., 1".18) in X and 
Y directions, and constructed light curves with photons 
from the sub-pixels (as in §2.1) within a 4" circular re¬ 
gion centered at each grid point. L-S periodograms were 
calculated as in §2.2, and the power for the strongest 
peak was chosen as the value for each grid point to con¬ 
struct the periodicity map as displayed in Figure [31(a). 

This periodicity map exhibits two strong periods at 
54.6 ±2.1 days and 62.0 ± 3.3 days (denoted by the red 
diamond and the blue triangle, respectively, in Figure 
|3)a)), which are positionally coincident with X-3 and X- 
4 respectively. In contrast, there is almost no significant 
periodicity near X-1. The individual periodograms of X- 
1, X-2, X-3, and X-4 with source regions of 4" in radii 
are displayed in Figure |3Kb). The highest peak in the 
periodogram of X-1 is at ~ 62.0 days, but with much 
lower confidence levels of 56% and 69% for white noise 


and red noise background respectively. Except for X- 
1, the periodograms of X-2, X-3, and X-4 all show the 
~ 55— day and ^ 62— day period, and the significance 
of the highest peaks in them are larger than 99% based 
on both white noise and red noise background. This 
implies that the ^ 55— and ^ 62— days peaks are gener¬ 
ated by real periodic signals with high probability rather 
than white or red noises, and the photons from the two 
signals are mixed up in the region containing X-2, X-3, 
and X-4 (dark part in Figure |3)a)), while they still can 
be distinguished from the counts of X-1. Thus the 62- 
day periodicity is most likely associated with the region 
4" southeast of X-1, instead of X-1 as expected. 

3. DISCUSSION 

By combining the RXTE and Swift light curves, we 
analyzed a longer term light curve. The analysis confirms 
phase changes of the 62-day period of M82, indicating 
that this period is not stable. A detailed timing analysis 
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Fig. 3.— (a) Periodicity map showing the L-S periodogram maximum power for a grid with steps of 0.5 pixels (i.e., 1".18). The 

periodograms were computed with photons from subpixels within a 4” circular region centered at each grid point. The red circle is centered 
on X-1 with a radius of 4”. The red diamond and blue triangle mark the center of the approximate 55- and 62- day periodicities respectively. 
It is clear that there is no significant periodicity within the 4"circular region of X-1. (b) L-S periodograms of X-1 (black filled line), X-2 
(green line), X-3 (red line), and X-4 (blue line) extracted from source region of 4" in radii. It shows no significant periodicity in X-1, 
whereas X-2 and X-3 show two significant periodicities at ~ 55 and ~ 62 days, and X-4 shows the strongest period at 62 days. The 
confidence levels of the highest peaks in X-1, X-2, X-3, and X-4 based on red noise are 69%, 99.6%, 99.7%, and 99.8% respectively. Dashed, 
horizontal lines have the same meaning as in Figure2 (a). 


of M82 shows two significant periodicities at 55 and 62 
days in a region 4" southeast of X-1, which contains X-2, 
X-3, and X-4. Given the X-ray localization precision of 
XRT, the periodicity map can hardly resolve X-2, X-3, 
and X-4 clearly, but it can easily distinguish the three 
ULXs from X-1. This indicates that the 62-day period is 
not from M82 X-1, but most likely from the three ULXs. 
The instability of the 62-day period may be due to the 
comprehensive effect of several different periodic X-ray 
sources, instead of a precessing accretion disk. 

It is not yet clear which sources produce the ^55— 
and ~ 62— day periodicities. However, X-4 can be ex¬ 


cluded first since it is a luminous young X-ray supernova 
remn ant candidate (iKong et al.l 120071 iChiang fc Kond 
[MU) , whose flux cannot fluctuate periodically. X-2 
(also referred to as X42.3-1-59) is interpreted as a mag¬ 
netized n eutron star system w ith an orbital period of 
2.5 days (iBachetti et al.ll20T^ . Superorbital periods of 
neutr on star binaries can ra nge from 25 days to 160 
days (jWiiers fc Pringle! 119990 . The ~ 55 days may 
be a superorbital period of X-2, but it may be not 
yalid because the X-ray luminosity of X-2 is seyeral 
orders of ma gnitude larger than f or the other neutron 
star binaries (iMa tsuba et al.lll995: [Reynolds et al.lll997t 


lOrlandini et al.lll998HSuchv et al.ll2OO80. 

For X-3, it is suggested as a background actiye galac- 
tic nucleus becau se of its soft excess below 2 keV 
(|Kong et al.l I2OO70 . Howeyer, soft excess is also com¬ 
mon in an obscured high-mass X-ray binary (HMXB; 


e.g., IGR J19140-F0951 and IGR .116207-51 29: IPrat et a.1.1 
120081 iTomsick et al.l 120090. X-3 h as also been regarded 
as an HMXB ( Mineo et al.l I2O120 . Given the strong 
star-forming activities in M82, HMX Bs are expected 
to do minate the bright X-ray sources (|Chiang fc Kongl 
ImH). Assuming the isotropic unabsorbed emission, the 
maximu m 0.5-10 keV lumi nosity of X-3 is 1.7 x 10^° 
erg s' ^ (|Kong et al.l I2OO70 . which is approximate to 
the (maximum-minimum) luminosity modulation of M82 
(^ 1.6 X 10^° erg s“^ ; KF07). Therefore, X-3 is bright 
enough to dominant the periodicity in M82 if X-3 is an 
HMXB with an orbital period of 62 days. 

In summary, the limited spatial resolution of XRT is 
the bottleneck in this study, but it is already indicating 
that the 62-day periodicity of M82 may not come from 
X-1. To obtain a reliable location of the periodicities, 
more monitoring of M82 is needed. 
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